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Abstract 

The  initial  vibrational  energy  distribution  of  molecules  that  collide 
with  a  hot  surface  has  been  varied.  The  effect  of  the  variation  on  the 
collisional  reaction  probability  (Pj  has  been  studied  under  single  collision 
conditions.  These  experiments  provide  a  more  sensitive  test  of  the  relative 
suitability  of  various  analytical  forms  for  the  collisional  transition  proba¬ 
bility  matrix,  P  ,  than  has  been  achieved  previously  in  earlier  VEM  experi- 
ments.  The  reaction  system  is  the  isomerization  of  cyclobutene  to  1,3  buta¬ 
diene.  A  seasoned  fused  quartz  surface  was  used  over  the  temperature  range  Tr  = 
600  K  -  900  K.  Variation  of  the  initial  vibrational  energy  population  vector  of 
cyclobutene  molecules  was  made  by  change  of  their  initial  temperature  Tc  in  the 
range  273  K  -  620  K.  The  experimental  collisional  efficiency,  81  .  declined 
from  0.31  to  0.013  over  the  combination  temperature  ranges  Tr,Tc  »  600,500 
to  900,273.  Stochastic  models  of  the  vibrational  transition  probability  were 
fitted  to  the  data  and  provided  values  for  the  average  amount  of  energy 
(<AE'>r  )  transferred  from  the  hot  molecules  in  a  down  transition  from  the 
threshold  energy  level  EQ;  Gaussian  or  Boltzmann  Exponential  forms  of  P  prove 

to  be  the  most  suitable  to  fit  the  data.  Calculated  values,  <AE ' >F  on  the 

-1  0  -1 

basis  of  a  Gaussian  function  model  for  P  declined  from  5600  cm  to  3700  cm 

* 

with  increase  in  the  surface  temperature  from  600  K  to  900  K. 
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INTRODUCTION 

Because  of  the  Interest  in  gas-surface  interaction*  energy  accommodation 
at  a  surface  continues  to  be  actively  studied  both  experimentally  and  theo¬ 
retically.^"5  Very  little  work  has  been  done  on  colllsional  transfer  of 
vibrational  energy  between  a  surface  and  polyatomic  molecules  at  the  level  of 
vibrational  excitation  corresponding  to  reaction  in  non-catalytic  processes. 

A  technique  for  this  purpose  called  the  Variable  Encounter  Method  (VEM),^  has 
been  reported  recently.  VEM  provides  a  measure  of  vibrational  energy  trans¬ 
fer  between  initially  cold  molecules  and  a  hot  surface:  gas  molecules  describe 
a  random  walk  along  an  energy  coordinate  until  they  reach  an  absorbing  level  - 
the  critical  reaction  threshold  EQ  characteristic  of  the  homogeneous  gas 
reaction.  Such  relaxation  of  non-equilibrium  vibrational  energy  distributions 
was  described  theoretically  a  number  of  years  ago  by  Rubin  and  Shuler,6  Kim,7 

O 

and  by  Widom. 

4  9-12 

In  the  previous  VEM  work,  ’  gas  molecules  were  thermally  equilibrated 
in  a  reservoir  flask  at  some  low  temperature  and  then  entered  a  hot  reactor  in 
which  they  experienced  a  known  (and  experimentally  variable)  average  number, 
m,  of  sequential  collisions  with  the  surface  before  leaving  the  reactor  and 
re-equilibrating  to  their  initial  low  temperature;  m  was  varied  between  2  and 
30.  A  fraction  of  the  molecules  are  excited  above  the  reaction  threshold 
energy  during  the  series  of  wall  collisions  and  experience  homogeneous  uni- 
molecular  reaction.  The  nature  of  the  collislonal  relaxation  of  the  vibrational 
energy  transients  may  be  deduced  from  the  experimental  reaction  rate.  Syscems 
involving  surface  catalysis  have  been  avoided. 

Recently,  single  collision  measurements  (m-1)  by  the  VEM  technique  have 
been  described.  In  the  present  paper  we  iia'*e  applied  the  single  collision 
condition  to  the  study  of  the  effect  of  variation  of  the  initial  vibrational 
energy  on  the  single  collision  reaction  (transition)  probability v> The  surface 
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Is  a  seasoned  fused  quartz  finger.  Although  such  a  surface  Is  not  well 
defined,  It  Is  the  conventional  experimental  surface  of  thermal  unlmolecular 
kinetics.  Unlike  low-energy  rion-reactlve  accommodation  studies  that  Involve 
translational,  rotational  and  vibrational  degrees  of  freedom,  vibrational 
energy  Is  the  main  contributor  to  the  reaction  probability  for  reacting  sys¬ 
tems  that  Involve  small  or  zero  centrifugal  energy  barrier  effects;  In  such 

14 

case,  rotational  energy  plays  little  role.  This  Is  the  case  here  for  the 
study  of  cyclobutene  unlmolecular  Isomerization  to  1,3-butadiene.  The  reaction 
has  a  low  activation  energy,  32.4  kcal  mole”1,  and  proceeds  with  no  evidence 
for  surface  catalysis.10,11,13 

The  vibrational  energy  population  vector  of  molecules  after  one  collision, 
N,  ,  is  given  by  N,  =  PN  ,  where  P  Is  a  colllsional  transition  probability 
matrix,  and  N is  the  Initial  vibrational  energy  population  vector  that 
corresponds  to  the  thermal  Boltzmann  distribution  at  the  low  temperature  of  the 
gas  reservoir  wall.  Since  no  comprehensive  theory  of  gas-surface  colllsional 
interaction  exists,  at  le2st  for  complex  molecules  on  these  surfaces,  down- jump 
transition  probabilities  were  constructed  according  to  various  assumed  trial 
analytical  models.  Four  models  have  been  applied  to  the  present  interpre¬ 
tations;  exponential  (E),  Gaussian  (G),  and  Boltzmann  distribution  weighted 
exponential  (BE)  and  Gaussian  (BG)  functions  have  been  used  in  conventional 
thermal15  and  VEM10,13  systems.  Up-transition  probabilities  were  constructed 
from  the  down-jump  transitions,  with  use  of  detailed  balance  and  completeness. 
Vibrational  energy  transfer  above  the  reaction  threshold,  E.  +  E,  (E.  £  E  ), 
is  the  measured  process  of  the  reaction  in  a  single  collision  reactive  system. 
The  transition  probability  matrix  Pisa  function  of  the  reactor  finger 
surface  temperature.  For  simplicity,  it  was  taken  to  be  independent  of  the 
temperature  of  the  reservoir  wall  and,  hence,  of  variation  of  mean  translational 
energy  of  the  molecules  and,  therefore,  of  small  variation  of  the  residence 


time  of  colliding  molecules  on  the  reactor  surface. 

Variation  of  the  Initial  population  vector  by  change  of  the  wall  tempera 
ture  constrains  the  models  of  the  collisions!  probability  matrix  that  can 
reproduce  the  experimental  reaction  probability  function  and  provides  a  more 
sensitive  test  of  the  aptness  of  the  various  trial  functions  than  has  been 
available  hitherto  by  this  technique. 
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EXPERIMENTAL 


The  cyclobutene  {Columbia  Chemicals,  >  99.9%  )  showed  no  Impurities 
on  gas  chromatographic  analysis.  The  reaction  system  apparatus  used  here 
was  similar  to  that  of  the  previous  single  collision  experiments  (Flg.l  of  ref  13)  . 
The  reaction  vessel  consisted  of  a  3-t  spherical  pyrex  reservoir  flask  that 
was  provided  with  an  Internally  heated,  central  fused  quartz  finger.  The 
flask  and  the  quartz  finger  were  heated  Independently.  The  heater  control  of 
the  finger  Is  described  In  the  previous  work  .  The  wall  temperature  of 

the  flask  was  controlled  by  an  air  oven  consisting  of  a  30-cm  cubic  insulated 
box  equipped  with  two  plate  heaters  and  an  air-mixing  stirrer;  except  that  the 
lowest  flask  temperature  (273  K)  was  attained  by  cooling  the  flask  in  an  ice 
bath.  Temperature  measurements  were  made  with  four  chromel-alumel  thermo¬ 
couples  cemented  to  the  outside  of  the  reservoir  flask.  The  surface  tempera¬ 
ture  of  the  reactor  finger  was  varied  from  600  K  to  900  K,  while  the  wall 
temperature  of  the  flask  reservoir  was  varied  from  273  K  to  600  K.  Temperature 
deviation  of  the  finger  surface  was  a  maximum  of  ±  5°K  at  900  K,  and  that  of 
the  flask  wall  was  ±  10°K  at  600  K.  Before  kinetic  measurements  were  made, 
the  reactive  surface  of  fused  quartz  was  "seasoned"  at  the  highest  temperature 

and  the  seasoning  was  maintained  by  exposure  for  a  few  hours  to  cyclobutene  at 

-3 

a  pressure  of  %  5x10  torr  prior  to  each  run. 

The  system  was  run  in  both  static  and  flow  modes  in  the  pressure  region 
-4  -3  -4 

between  10  and  10  torr,  usually  2x10  torr.  for  the  static  mode,  the 
reactor  was  pumped  to  <  10~^  torr  before  substrate  was  Introduced.  The  reaction 
gas  was  transferred  for  analysis  after  a  predetermined  run  time.  For  the  flow 
mode,  the  Impedance  of  a  variable  exit  valve  was  calibrated  by  pressure-drop 
measurements.  Typical  residence  times  were  10  to  40  sec.  The  run  was 
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terminated  by  closing  the  variable  flow  valve  after  about  30  to  100  sec 

of  reaction  time.  More  detailed  description  of  the  flow  analysis  was  given 
)  <» 

earlier.  Pressure  measurements  were  made  with  an  MKS  Model  170  capacitance 
manometer. 

Product  analysis  was  performed  by  gas-liquid  chromatography  on  a 
5  ft  x  3/16  In.  squalane  column  on  Chromosorb  P  at  0°C  with  FID  detection. 


/ 
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RESULTS 

Observed  rate  constants  for  cyclobutene  Isomerization  to  1,3-butadiene 
were  calculated  from  the  product  yield.  The  total  observed  rate  constant 
kt  was  a  sum  of  reaction  due  to  heating  at  the  retctor  finger,  kr  ,  and  at  the 
reservoir  wall,  k^  ,  given  by 

kt(Tf,Tc)  -  kr(Tr,Tc)  +  kc(Tc)  , 

where  T  and  T  are  the  temperature  of  the  finger  surface  and  of  the  reservoir 
r  c 

wall,  respectively,  k  proves  negligible  for  small  T  values.  Independent 

w  V 

measurements  of  the  k  (T  )  were  made  under  the  experimental  condition,  T  *  T  , 

o  c  r  c 

with  use  of  a  minor  area  correction  for  the  relative  surface  area  of  the 

reactor  finger  and  the  reservoir  wall  (1:13).  The  reaction  rate  constants 

k.  and  k  are  shown  in  Fig.  1.  The  rate  of  k  rises  rapidly  with  an  increase 
t  c  c 

in  Tc  because  of  the  larger  area  of  the  reservoir  wall.  Experimental  measure¬ 
ments  of  k  were  rejected  in  the  cases  when  correction  for  k  amounted  to  >  50% 
r  c 

of  the  total  rate  k^;  these  rejected  values  are  marked  with  an  X  in  Fig.  1. 

The  reaction  probability  per  collision,  P  (T  ,T  ),  was  calculated  from  the 

U  I  V# 

finger  surface  reaction  constant  kr (Tr ,l'c ) .  Relative  values  of  the  reaction 
probability  are  plotted  as  a  function  of  the  T  temperature  in  Fig.  2  on  the 

w 

basis  of  Pc (Tr , 273 )  =  1.  Error  bars  represent  experimental  uncertainty.  As 
is  reasonable,  increase  of  T  has  greater  effect  when  T  is  smaller.  Experi- 

L  i 

mental  values  of  the  reaction  probability  per  collision  are  listed  in  Table  1. 

The  experimental  uncertainties  in  P  are  estimated  to  vary  from  %  10%  at  low  T 

temperature  to  as  much  as  20%  at  high  temperature.  The  collisional  efficiency, 

,  defined  as  S-j  5  Pc(Tr,Tc)/P*c(Tr)  also  9iven  Table  1.  The  elements 

of  P*C(T  )  are  simply  the  Boltzmann  distribution  population  vector  for  a  uni- 
c  r 

molecular  reaction  In  the  low  pressure  regime.  They  were  calculated  with  use  of 

16  * 

the  vibrational  frequency  assignments  of  Elliott  and  Frey. 
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DISCUSSIONS 

The  following  experimental  and  theoretical  simplifications  In  the  present 
type  of  study  which  Increase  the  power  of  the  method  should  be  noted.  First, 
no  knowledge  of  molecular  collision  cross-sections  is  Involved.  Second,  If 
the  reaction  process  Is  complete  within  the  free  flight  time  of  the  molecule 
when  It  leaves  the  surface,  then  the  magnitude  of  the  specific  decomposition 
probabilities  k(E)  do  not  enter.  The  reaction  Is  in  the  unimolecular  low  pressure 
regime  and  neither  a  postulated  activated  complex  structure  nor  RRKM  calculations 
are  required.  This  Is  the  case  here;  decomposition  under  all  conditions  exceeds 
99%  of  the  activated  molecules. 

Variation  of  Pc(Tr,Tc)  with  Tc  Is  shown  in  Fig.  3  for  each  value  of  Tr 
studied.  Several  fits  to  the  experimental  data  were  made  on  the  basis  of 
various  assumed  analytical  forms.  The  mathematical  forms  are  exponential,  E  , 
gaussian,  G  ,  Boltzmann  exponential,  BE,  and  Boltzmann  gaussian,  BG.  The  corres¬ 
ponding  forms  of  the  probability,  P(AE),  of  a  down  transition  of  energy,  AE  , 
are  given  by  Eqs.  (1)  -  (4). 


t»E (AE)  • 

c^exp(-AE/<AE>) 

(1) 

PG(4E)  - 

c2exp(-(AE-AEmp)2/2a2) 

(2) 

>“(AE) 

»  c3B.PE(AE) 

(3) 

P^IAE) 

*  c4b.pg(ae) 

(4) 

Here,  <AE>,  AE  ,  and  a  are  adjustable  parameters,  the  c4's  are  constants, 
mp  • 

and  B^  =  g^exp(*E|/RTr)  is  a  normalized  Boltzmann  distribution  characteristic 
of  the  finger  surface  temperature,  Tr,  where  g^  Is  the  density  of  vibrational 
states  at  energy  level  E^ .  The  E  and  G  forms  of  Eqs.  (1)  and  (2)  have  been 
termed  "flat"  forms  previously  since  the  transition  probabilities  are  Indepen¬ 
dent  of  energy  level.  They  are  inadequate  to  describe  near-strong  collider 


\ 


behavior  since,  as  strong  collider  behavior  Is  approached,  the  transition 
probability  should  approach  the  Boltzmann  probability  function,  ^(AE)  ■ 
Hence,  the  Bo1t2mann  weighted  models  of  Eqs.  (3)  and  (4)  were  also  applied 
to  the  data.  In  earlier  VEM  studies  (m  >  1),  the  E  model  produced  the  best 
agreement  with  the  experimental  data  of  cyclobutane  over  the  range  of  reaction 

o 

surface  temperatures  from  749  K  -  1126  K  and  with  the  methyl  cyclopropane 

12 

system  from  800  K  -  1130  K.  The  G  model  was  found  to  give  the  best  overall 

4  4 

fit  to  the  data  for  1 .l-cyclopropane-d^i  while  the  BE  model  was  used  for  the 
cyclobutene  reactive  system  (m«l)  over  the  lower  temperature  range13  400  K  - 
915  K,  although  the  E  and  G  models  fit  better  at  higher  temperatures.11 

Results  of  least  squares  fitting  to  the  experimental  P  curves  for  each 
of  the  four  analytical  formsof  P  tested  are  also  shown  In  Fig.  3;  variation  of 
<AE>,  or  AEmp,  was  the  parameter  of  fit.  The  best  fit  values  are  listed  in 
Table  2.  We  use  and  exhibit  these  average  best-fit  values  rather  than  enumer¬ 
ating  a  multitude  of  individual  best-fit  values  for  each  T  for  a  given  T  . 

For  the  G  and  BG  models,  calculation  was  made  with  the  standard 


deviation  set  at  o  *  0.7  AEmp,  as  in  our  earlier  calculations.  The  calculated 

values  of  <AE‘>F  ,  the  average  amount  of  energy  transferred  in  a  down  transi- 
0 

tlon  from  the  threshold  energy  level,  EQ  ,  are  also  given.  Other  quantities 
of  Interest  are  AEay,  AE*v.  AEflV  is  defined  as  AEflV  =  Ef-§c  ,  where  Ec  is 
the  average  energy  of  molecules  at  temperature  T  before  collision  and  Ef  is 
the  average  energy  after  collision.  The  other  useful  quantity  AE*v  is  the 
amount  of  average  energy  transferred  above  the  critical  threshold  defined  as 


Mav  =  ^Ef^Ef  >E0‘  V* 
calculations  in  Table  3. 


They  are  listed  with  results  of  G  and  BE  model 


The  conventional  vibrational  energy  accommodation  coefficient,  a  ,  is 
defined  as  a  *  AEau/(E  -Er),  where  E„  is  the  average  vibrational  energy  of 
molecules  equilibrated  at  reactor  temperature  T^;  values  are  given  in  Table  3. 
Average  vibrational  energies  of  cyclobutene  molecules  for  the  varlpu,  thermal 
equilibrium  populations  are  given  in  Table  4. 
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The  features  of  the  experimental  P  curve  in  Fig.  3  are  relatively  flat 

c 

lines  at  Tc  <  400  K  and  a  steeper  rise  at  Tc  >  400  K.  The  P£  curves  calculated 
with  both  the  G  and  BE  models  agree  substantially  with  the  data,  with  the 
former  giving  somewhat  better  fit.  But  the  only  model  that  can  fit  the  In* 
efficient  flat  behavior  In  P  at  low  T  Is  the  E  model.  In  general,  however, 

U  V 

both  E  and  BG  curves  differ  most  from  the  bulk  of  the  data  and  Illustrate  the 
extreme  behavior  of  weak  and  strong  collider  behavior,  respectively. 

If  one  compares  AEau  for  G  and  BE  models  which  are  both  fairly  well  fitted 

iV 

to  the  data,  the  relative  values  seem  paradoxical:  that  for  G  Is  very  big  and 
that  for  BE  very  small;  although  both  reproduce  the  experimental  values  of  P 

v* 

and  the  general  shape  in  Fig.  3.  Again  a  is  large  for  G  and  small  for  BE 

(Table  3).  The  explanation  is  to  be  found  In  Fig.  4,  a  typical  plot.  The  plot 

G  G 

of  N,  follows  Ny  at  low  energies  corresponding  to  high  a  values.  However, 

r  r  BE 

above  EQ  the  Np  values  drop  relative  to  Nip  which  never  attain  Ny  but  do 
not  tail  off  as  fast  as  N,^  at  higher  energies.  Thus,  AE*  In  Table  2  is 

*•«  flV 

DC 

larger  for  BE  than  for  the  other  models.  So  although  a  is  small,  the  model 
gives  the  same  P  result  as  the  G  model  because  the  reaction  probability  is 
determined  by  the  fraction  of  the  original  population  transported  above  energy,  Eq 
The  P  matrix  may  be  partitioned  into  four  parts  at  the  energy  Zq  as  shown 
in  Eq.  (5).  For  the  consideration  of  energy  accommodation  coefficient  a,  the 


(5) 


transition  probability  involving  states  excited  below  EQ  in  the  sub-matrix 
Pj  provides  the  main  contribution  to  the  vector  leading  to  the  average 
energy  change  of  molecules,  AE.  The  population  of  reacting  states  is  derived 

dV 

from  the  transition  probability  to  states  above  Eq,  l.e.,  Pjj,  which 
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Is  connected  with  down- transitions  P^j  by  detailed  balance.  In  the  VEM 
method  with  m  =  1 ,  information  is  gained  regarding  elements  of  Pjj 
and  Pjjj  .  Little  information  is  obtained  about  region  Pj,  so  that 
the  models  that  fit  P  are  Inadequate  to  predict  a  values.  This  situation 

L 

Is  modified  for  m  >  1  where  more  Information  can  be  extracted  concerning 

region  I,  but  with  less  ability  to  discriminate  between  the  models  that  fit 

region  II.  Obviously,  a  and  P  measurements  are  complementary  for  the  elucida- 

tion  of  the  most  appropriate  analytical  forms.  Draper  and  Rosenblatt  measured 

vibrational  energy  accommodation  coefficients  of  hydrocarbon  on  a  "dirty"  metal 

2 

surface  at  room  temperature  by  a  vibrating  diaphragm  technique.  Although  those 

experiments  are  not  directly  comparable  to  ours,  their  values  of  ayib  are  0.7 

to  0.9  for  normal  paraffins  (C^-Cg).  Since  the  G  model  provides  calculated 

a  values  nearly  close  to  unity,  this  might  be  considered  to  be  the  preferred 

form  here.  However,  it  is  evident  that  no  single  analytical  form  so  far  tested 

gives  optimum  fit  of  the  data  over  the  whole  range  of  experimental  conditions. 

2 

The  recent  measurements  of  a  by  Foner  and  Hudson  for  butane  on  Pt  at  1000°C 
support  the  BE  model  instead. 

An  important  difference  exists  between  weak  collider  and  strong  collider 
behavior  in  a  m  =  1  system.  For  a  strong  collider,  activation  to  the  levels 
above  En  occurs  mainly  from  the  bulk  of  the  thermal  distribution  at  T  .  Only 
for  weak  colliders  does  the  contribution  come  mainly  from  the  sparsely  populated 
levels  just  below  Eq.  This  is  illustrated  in  Fig.  5  for  colliders  of  different 
strengths. 

The  decline  in  collisional  efficiency  that  occurs  with  rise  of  T  ,  at 

constant  T  in  (Table  1)  has  been  observed  in  the  previous  cyclobutene  VEM 
c 

studies.  The  gas-surface  interaction  approaches  strong-collider  behavior 

at  T  <  450  K.  This  trend  appears  to  be  related  to  the  fact  that  as  T  is 
r  r 

lowered,  the  time  that  the  molecule  spends  trapped  in  the  gas-surface 

It  13 

potential  well  Increases.  ‘ 
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Table  1.  Experimental  reaction  probability  and  colllslonal  efficiency 
In  the  m=l  cyclobutene  system. 


r\c'*' 

273 

400 

500 

600 

900 

P  xlO4 
c 

1.48 

1.54 

2.50 

6.50 

‘ia 

0.013 

0.013 

0.022 

0.056 

800 

Pcxl05 

4.18 

4.40 

8.50 

24.8 

*1 

0.019 

0.020 

0.039 

0.113 

700 

P  xlO6 
c 

9.42 

10.8 

20.5 

el 

0.044 

0.051 

0.096 

600 

P  xlO7 
c 

9.33 

11.3 

23.5 

S1 

0.123 

0.149 

0.310 

colli sional 

efficiency  ^  is 

defined  as 

61  ■  Pc<Tr 

,Tc>/p 

sc, 


where  P^C(T  )  is  the  reaction  probability  due  to  a  strong  collider. 

v  r 

It  is  the  analog  of  the  conventional  homogeneous  unimolecular 
collision  efficiency,  6  . 
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Table  2.  Energy  transfer  parameters  (cm’1)  for  various 
models  in  the  m  = 1  cyclobutene  system. 


Tr 

Model 

<AE>  or  AEmna 
mp 

<a£,>Ed 

(<>U 

900 

2450 

2270d 

840 

800 

3260 

2780d 

770 

700 

E 

4080 

3 1 20  d 

660 

600 

b 

— 

— 

900 

3380 

n 

3740  w 

420 

800 

3700 

4090  6 

440 

700 

G 

4200 

4590® 

450 

600 

5630 

5640 

480 

900 

1470 

2850 

1110 

800 

BE 

1340 

3490 

950 

700 

1200 

4450 

780 

600 

1150 

6290 

610 

900 

2880 

4390 

380 

800 

8G 

2850 

4900 

360 

700 

3110 

6070 

400 

600 

3290 

7300 

410 

a) 

For  G  and  BG  models,  the  calculations  were  made 

with  the  standard 

deviation  set 

at  o  =  0.7  AE  . 

mp 

b) 

Least  square 

fit  was  impossible  for 

any  physically  sensible  <AE>  value. 

c)  (AE^),.  a  (E*-En)r  *  values  given  for  case  T  =  400  K. 

aV  T  U  t^tg  C 

d)  Differ  from  column  3  due  to  truncation  effects  approaching  E = 0  where 

i' 

the  transition  probabilities  were  constrained  not  to  exceed  Boltzmann. 

e)  Differ  from  column  3  due  to  truncation  near  E*Eq  as  well  as  E*0.  (See  (d )) 


Table  3.  Average  energy  transfer,  AEflV  (cnf^),and  vibrational 
accommodation  coefficients,^,  in  the  m  *  1  cyclobutene  system. 


TC(K) 


Tr  (K) 

Model 

273 

400 

500 

600 

900 

3850 

0.87 

3400 ( 420 )b 

0.86 

2880 

0.84 

2230 

o 

*  i 

OO  i 

.'•'0 

800 

3200 

0.94 

2730(440) 

0.92 

2170 

0.90 

1500 

0.87 

700 

b 

2420 

0.97 

1940(450) 

0.95 

1380 

0.92 

720 

0.89 

600 

1620 

0.96 

1150(480) 

0.93 

620 

0.90 

900 

770 

0.17 

740(1110) 

0.19 

690 

0.20 

590 

0.21 

800 

BE 

690 

0.20 

650(950) 

0.22 

560 

0.23 

430 

0.25 

700 

610 

0.24 

530(780) 

0.26 

420 

0.28 

240 

0.30 

600 

560 

0.33 

430(610) 

0.35 

260 

0.37 

a)  AEau 

=  'Ef 

Ec;  a 

=  AEav/(E 

-s 

1 

mi 

o 

b)  Values  of  (<v)Ef>EQ 


in  parenthesis  from  Table  2  shown  for  comparison. 
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Figure  Captions 

Fig.  1  Observed  rate  constants  for  cyclobutene  Isomerization  to 
1,3  butadiene:  the  open  symbols  are  for  the  total 
activation  reaction,  l^{Tr,Tc),  and  the  closed  circles  are  for 
the  reservoir  wall  reaction,  k  (T  ).  The  X's  represent  data 

V  V 

which  was  dropped  for  reason  of  the  large  correction  on  k(T  ,T  ) 
due. to  kc(Tc). 

Fig.  2  Relative  reaction  probability  per  collision  on  the  basis  P  (T  ,273)  *  1. 

V  I 

Fig.  3  Results  of  the  least  square  fit  to  the  experimental  P  curve  for  each 

finger  surface  temperature:  experimental  curve  ( - );  E  ( - ); 

G  ( - );  BE  ( - );  and  BG  ( - ).  The  E  curve  could  not  fit 

the  data  for  T  *  600  K  and  Is  not  shown, 
r 

Fig.  4  Vibrational  energy  population  of  cyclobutene  for  the  reaction  con¬ 
dition  (900,  400).  N®q  and  N*q  are  the  thermally  equilibrated 

F  w 

Boltzmann  populations  at  the  temperatures  T  and  T  ,  respectively. 

*  V 

k(E)  is  the  isomerization  rate  constant  calculated  by  RRKM  theory. 

BE  and  G  represent  the  population  after  one  collision, calculated 
from  the  least  square  fit  with  the  BE  and  G  models.  These  normalized 
population  curves  refer  to  the  two  left  ordinates. 

Fig.  5  The  relative  probability  of  up-transitions  to  the  population  at 

En,  given  by  pP  .•  n  .  ,  is  illustrated  for  G  and  BE  models  with 

two  sets  of  values  of  the  parameters,  corresponding  to  strong(er)  and 

weak(er)  collisions.  The  several  sets  for  the  two  models  give  roughly 

similar  values  of  P  .  The  calculation  was  made  for  the  condition 

c 

T  *  900  K  and  T  -  400  Kj  <AE>(cm_1)(P  )  =  4000  (8.53x1 0-4)  and 

i  C 


■-  Vv  18 

1 500(3.44x1 0“9)  for  the  G  model;  and  2000  (7.12xl0“4)  and 
500  (2.94xl0‘9)  for  the  BE  model.  The  ordinate  Is  In  arbitrary 
linear  units  rather  than  log  units  so  as  to  bring  out  detail. 
Some  Irregularity  due  to  quantized  nature  of  the  state  density 
Is  revealed  at  E  <  2000  cm"1. 
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